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Nonlocal Manipulation of Dimer Motion at Ge(OOl) Clean Surface via Hot Carriers 

in the Surface States 

Yasumasa TAKAGI*, Yoshihide YOSHIMOTO, Kan NAKATSUJI and Fumio KOMORJt 

Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa-shi, Chiba 277-8581, Japan 

Nonlocal one- dimensional motions of a topological defect are induced by electron tunneling 
through the dangling-bond states on the clean Ge(OOl) surface using scanning tunneling mi- 
croscopy below 80 K. The direction of the motion depends both on the energy of the carriers 
in the surface state and on the distance between the defect and the tunneling point. The re- 
sults are interpreted using an electronic excitation model by hot carriers injected to the surface 
states. The critical distance of the motion is anisotropic and consistent with the band structure 
of the surface states. 
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Local electron tunneling often induces motions of 
atoms and molecules at surfaces, and their dynamics 
has been studied using scanning tunneling microscopy 
(STM). In most of the cases, the electron is directly 
injected to or extracted from their localized electronic 
states just under the tip, and the electronic energy is 
transferred to the kinetic energy of atoms or molecules by 
inelastic tunneling processes. 1 3 On semiconductor sur- 
faces, the lifetime of the surface electron is longer than 
on metal surfaces, and a nonlocal excitation through a 
surface state has been expected. Previously, nonlocal mo- 
tions of the halogen adsorbates by the tunneling were re- 
ported on Si(lll) 4 and Si(OOl) 5 surfaces. However, the 
propagation of the injected hot carrier has not been stud- 
ied in detail. It is important to understand the mecha- 
nism of such nonlocal manipulations both for fundamen- 
tal physics on electron and lattice dynamics and for their 
application to nanostructure fabrications. 

In the present letter, we show motions of a topological 
defect on the clean Ge(001)-c(4x2) and -p(2x2) surfaces 
at 80 K as functions of the distance and the direction be- 
tween the tunneling point and the defect (cf. Fig. 1). The 
defect can be easily created and annihilated by tunneling 
current, 6 and thus is suitable for the quantitative study 
of the nonlocal manipulation. Its motion is anisotropi- 
cally induced by electron injection to the surface while 
the excitation is isotropic in the case of hole injection. 
The critical distance of the motion can be more than 
80 nm by the electron injection. Moreover the direction 
of the motion depends on the sample bias voltage and 
the local superstructure. The observed anisotropy of the 
defect motion can be understood by considering the dis- 
persion relations of the surface electronic bands. 

The Ge(001) clean surface consists of the buckled 
dimers, and its buckling orientation is alternate in the di- 
rection perpendicular to the dimer axis. We call the row 
of the dimers in this direction a "dimer row" . There are 
two stable arrangements of the dimers; one is the c(4x2) 
structure, and the other is p(2x2). Figures l(a,b) show 
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Fig. 1. (color online) (a,b) STM images including topological de- 
fects, A, B, which were created on a c(4x2) surface (a), or on 
a p(2x2) surface (b). Protrusions at C in (a) and D in (b) are 
impurities as fixed markers, (c) A schematic model of the defect 
between the two superstructures, (d) STM image after the bias 
pulse from - 0.4 V to 0.8 V for 1 sec at P in (a). The STM tip 
position is fixed during the pulse. 



STM images at 80 K including a narrow p(2x 2) area on 
c(4x2) surface (a), and a c(4x2) area on p(2x2) surface 
(b). At the boundary between the two structures in the 
dimer row, there are two adjacent dimers with the same 
buckling orientation as a topological defect. Its model is 
schematically illustrated in Fig. 1(c). It is imaged higher 
than the other Ge atoms on the surface when the sample 
bias voltage Vb is negative. Here, we call it a "kink" . For 
simplicity, we also define the kink shown in Fig. 1(a) as 
A- type and that in Fig. 1(b) as B-type. 

The STM images were observed at 80 K in an ultra- 
high vacuum system. 7 The c(4x2) Ge(OOl) clean surface 
was obtained by several repetitions of Ar ion sputtering 
and annealing. We made a large area of the p(2x2) struc- 
tures by scanning the surface with V& = 1.2 V and the 
tunneling current I t = 1 nA, 6 ' 8 typically. We can con- 
vert the p(2x2) surface to c(4x2) again by scanning the 
surface with Vb = — 1.2 V. 

The kink is created and moved on the surface by inject- 
ing electrons or holes at its neighboring point using STM 
below 80 K. For example, the A- type kink was made by 
changing Vb from — 0.4 V to 0.8 V for 3 sec with fix- 
ing the tip position over a dimer in the c(4x2) surface. 
Generally, a pulse operation creates two kinks, and there 
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forms a one-dimensional area of p(2x2) between them. 
In the present study, the one end of the area was fixed at 
a step edge, and we focused on the motion of the kink on 
the terrace. The surface can be clearly imaged without 
moving the kink using Vb = — 0.4 V and I t = 1 nA. 
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Fig. 2. The distance between the kink and the tunneling point in 
the dimer-row direction after each bias pulse from — 0.4 V to 
0.8 V (a), 0.9 V (b) and 0.7 V (c) for 1 sec was applied at the 
point on the dimer row including the kink. The current during 
each pulse was 25 nA (a), 33 nA (b) and 17 nA (c). 



A Vb pulse at the position of P on the surface shown in 
Fig. 1(a) randomly moves the A- type kink in the dimer- 
row direction without any other changes in the surface 
structure as in Fig. 1(d). We observed STM images with 
Vb = — 0.4 V after each pulse to 0.8 V for 1 sec. Figure 
2(a) summarizes the random motion of the kink posi- 
tion by the pulse along the dimer row. The center of the 
random motion is 10 nm apart from the position of the 
tunneling. The separation between the center and P in- 
creases with increasing Vb during the pulse. It is 30 nm 
with the pulse to 0.9 V as in Fig. 2(b), and over 50 nm 
with the pulse to 1.0 V. On the other hand, when we 
use the pulses to 0.7 or 0.6 V, the kink first comes close 
to the tunneling point, and moves randomly around the 
point as in Fig. 2(c) for the 0.7 V pulse. The width of the 
random motion by the 0.6 V pulse is smaller than that 
by the 0.7 V pulse, and the positive pulse applied to Vb 
of less than 0.5 V never induces the kink motion. In the 
case of the B-type kink, the distance between the tun- 
neling point and the center of the random walk rapidly 
increases with increasing Vb of the pulse above 0.7 V, and 
becomes more than 50 nm with the pulse to 0.75 V for 
1 sec. 

Random motion of the kink by a negative bias pulse 
was previously reported. 6 The distance between the cen- 
ter of the kink motion and the tunneling point increases 
with decreasing Vb of the pulse below — 0.7 V. In the 
negative bias case, there is no bias range where the kink 
comes back to the tunneling point; changing Vb of the 
pulse to above — 0.6 V does not move the kink. The pre- 
ferred direction of the kink motion, thus, largely depends 
on the bias voltage and the distance from the tunneling 
point. 
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Fig. 3. (color online) (a,b) Probability of the kink motion per bias 
pulse as a function of the distance between the tunneling point 
for the A-type (a) and B-type (b) kinks. The pulse is from — 0.4 
V to 0.7 V for 1 sec. The dotted curves are obtained using eq. 1 
with the parameters fitted to the data, (c) Two-dimensional plot 
of the tunneling point where the probability of the kink motion 
is 0.5 with positive pulses from — 0.4 V to 0.6 and 0.7 V for 
the A-type kink, and to 0.7 V for the B-type kink. The kink 
is located at the origin, (d) The same plot with negative pulses 
from - 0.4 V to - 0.7, - 0.8, and - 0.9 V for the A-type kink. 



Motions of the kink were repeatedly observed as func- 
tions of the distance and direction from the tunneling 
point for the both types of the kink. We show in Figs. 
3(a,b) the probability of the motion per pulse vs. the dis- 
tance in the dimer-axis direction. The STM tip position 
is fixed with Vb = — 0.4 V and I t = 1 nA, and the pulse 
is to 0.7 V for 1 sec. In this condition, we injected 1.1 
x 10 11 electrons per pulse at the tunneling point. The 
probability of the kink motion rapidly decreases at a cer- 
tain distance. We define here the critical distance of the 
motion, L c , as the distance where the probability is 0.5. 
It depends on the type of the kink, the bias voltage, and 
the direction from the kink to the tunneling point. 

In Fig. 3(c), we two-dimensionally plot the points 
where the probability is 0.5 for the both types of the 
kink. The total number of injected electrons for the 0.6 
V pulse was 8 x 10 10 . The critical distance is highly 
anisotropic and exceeds 80 nm in the dimer-row direc- 
tion even for the pulse to 0.6 V. Thus, the ratio of L c 
in the dimer-row direction to that in the dimer-axis di- 
rection is more than ten. On the other hand, the ratio is 
two in the case of the hole injection to the surface as in 
Fig. 3(d). The total numbers of injected electrons were 
2.5 x 10 10 for the - 0.7 V pulse, 4.4 x 10 10 for - 0.8 V, 
and 6.3 x 10 10 for - 0.9 V. 

The critical distance observed in the dimer-axis direc- 
tion by electron injection depends on the type of the kink. 
To understand the origin of the difference, we prepared 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 3 



thin c(4x2) areas in different ways between the B-type 
kink and the tunneling point on the p(2x2) surface as 
in Figs. 4(a-d). We measured L c on these surfaces using 
the same pulse from — 0.4 to 0.7 V. Here we always in- 
jected the current on the p(2x2) surface to maintain the 
tunneling condition the same. We found little difference 
of L c between the surfaces shown in Figs. 4(b) and 4(c). 
These two surfaces have the same width of c(4x2) area in 
total, but the arrangement of the two areas is different. 
This means that L c depends mainly on the total width 
W. The observed L c is almost a linear function of W as 
in Fig. 4(e). 
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Fig. 4. (color online) (a-d) STM images showing the surfaces with 
different local structures used for the measurements of the IB- 
type kink motions. The total width of the c(4x2) area, W, is the 
same in (b) and (c). (e) The critical distance L c in the dimer-axis 
direction as a function of W . The slope of the dotted line is — 
1.5. 

For the discussion of the experimental results, we 
adopt the model that the kink motion is induced through 
an inelastic scattering process by the injected hot carri- 
ers in the substrate as in the case of the adsorbate mo- 
tions. 4,5 In the present case, the hot carrier maintains 
its energy enough to induce the kink motion during the 
propagation from the tunneling point to the kink in the 
surface state. The threshold energy of the kink motion 
is around 0.6 eV for both electron and hole injections. It 
is larger than the theoretically estimated barrier of the 
kink motion, 0.3 eV. 9 The apparent discrepancy, how- 
ever, can be removed if we consider the very small rate 
of the inelastic process through electron-phonon interac- 
tion below the observed threshold energy. 10 

The anisotropy of the critical distance depends on 
the electronic structures of the surface. The valence and 
conduction band structures of the Ge(001) surface have 



been studied by photoemission, inverse photoemission 
and tunneling spectroscopies. 11 14 The filled tt band is 
a nearly isotropic surface state, and is a resonance to 
the bulk valence band around T. On the other hand, the 
empty 7r* band is highly anisotropic, and is isolated from 
the bulk states in the band gap. It strongly disperses 
in the dimer-row direction, and there is little dispersion 
in the dimer-axis direction. These experimental results 
are qualitatively consistent with the band calculations 
for the c(4x2) surface 14 and for the 2x1 surface with 
the buckled dimers. 15 Consequently, the hole injected to 
the 7r band two-dimensionally propagates in the surface. 
This explains the observed nearly isotropic L c in Fig. 
3(d) while the anisotropic L c in Fig. 3(c) for the elec- 
tron injection is attributed to the anisotropy of the hot 
electron propagation in the 7r* band. 

The relaxation time of the hot carriers injected to the 
surface bands is crucial in the present model. It looses the 
energy through electron-phonon and electron-electron in- 
teraction. The dynamics of the hot electron in the 7r* 
band has not been studied for the clean Ge(001) sur- 
face, but was measured by time-resolved two-photon pho- 
toemission spectroscopy for the Si(001) surface, 16 which 
has a similar electronic structure to the Ge(001) surface 
around the Fermi energy. The hot electron in the 7r* band 
of the Si(001) surface inelastically scatters within the 
band via phonon emission to the band bottom at T in 1.5 
ps. After that, it further relaxes from the band bottom 
to the other states including the bulk states. Further- 
more, a time scale of 50 fs was reported for the intra- 
band inelastic scattering time via phonon emission using 
a five- wave- mixing method. 17 The short scattering time 
was explained by a strong electron-phonon interaction at 
the surface. On the Ge(001) surface, a similar intraband 
relaxation process of energy is expected to dominate the 
relaxation to the bulk by considering the resemblance of 
the 7r* band. In the tt band, the carrier relaxation time 
to the bulk states should be shorter than that in the 7r* 
band because the tt band is a resonance to the bulk states 
around T. 

When we assume the following semi-classical model of 
the hot electron propagation in the dimer-axis direction, 
we can reproduce the observed linear relation between 
L c and W . The validity of the model means that the 7r* 
electron has a quantum coherence length shorter than a 
few nm in the dimer-axis direction. The phonon emission 
with the short interval may destroy the coherence during 
the slow propagation of the hot electron in this direction. 

For the analysis, we simply suppose that the hot elec- 
tron loses energy to less than the threshold of the kink 
motion with a certain probability q during the electron 
transfer from a dimer row to the adjacent dimer row. Fur- 
thermore, we presume that the q for the c(4x2) structure, 
<2c(4x2) is different from that for p(2x2), q p (2x2)- Now, we 
consider the mixed surface structure which consists of 
m c ( 4X 2) dimer rows of the c(4x2) structure and m p ( 2X 2) 
dimer rows of the p(2x2) structure between the kink and 
the tunneling point as in Fig.4(a-d). Then, the probabil- 
ity P(m c ( 4x2 ), m p (2x2)) that at least one electron moves 
the kink through (ra c ( 4x2 ) + ^(2x2)) dimer rows after 
injecting N electrons can be written as 
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P( m c(4x2), m p(2x2)) 
= 1 - (1 - 5(1 - 4c(4x2)) mc(4x2) (1 - ^(2x2)) mp(2x2) ) 



(1) 



Here, s is the probability that a single hot electron at the 
kink moves it. We used the assumption that each motion 
of the hot electron between the dimer rows is independent 
in the semi-classical model. The dotted curves in Figs. 
3(a,b) are those fitted to this equation using the following 
parameters; s = 8.7 x 10 -9 , q c (4x2) — 0-42 and m p ( 2X 2) 
= for (a) and s = 2.2 x 10 -8 , m c ( 4X 2) = and q p (2x2) 
= 0.21 for (b) with the fixed TV = 1 x 10 11 for the both. 
It is noted that the parameter s should depend on the 
type of the kink. 

To show the linear relation between L c and W, we 
define m c = ra c ( 4x2 ) and m p = m p ( 2X 2) a t L C1 that is, 
P(m c , nip) — 0.5, L c = a(m c + m p ), W = am ci and 
a is the period of the dimer row. Then, the relation, 
^c/^co + rrip/mpo = 1, is obtained after calculations 
using eq. 1. Here, m c o and m p o are experimentally de- 
termined on the pure c(4x2) and p(2x2) surfaces by 
the equations, P(m c o,0) = 0.5 and P(0,m p o) = 0.5. 
The linear relation can be written as L c = — (m p o — 
m c o)W/m c o + am P Q. Using the observed values am c o = 
10.4 nm, and am p o — 27 nm in Figs. 3(a,b), we ob- 
tain, L c = — 1.6W + 27 nm. The coefficient of the linear 
term agrees with the slope of the dotted line in Fig. 4(e) 
while the constant term is higher than the experimental 
value. The discrepancy is attributed to the simple as- 
sumption in the semi-classical model. When the width of 
the c(4x2) area is small, a quantum coherence of elec- 
trons may not be neglected during the propagation in 
the dimer-axis direction. 

The reason for the preferred direction of the kink mo- 
tion is not clear. The results shown in Fig. 2 indicate that 
the direction depends on the distance from the tunnel- 
ing point and the bias voltage. The interaction between 
the electric dipole of the Ge dimer and the electric field 
due to Vb is a possible origin of the preference because 
it depends on the distance from the tunneling position. 
The kink motion is essentially the inversion of the buck- 
ling orientation of the Ge dimer, and the electric dipole 
energy in the electric field can exceed the energy differ- 
ence per dimer between p(2x2) and c(4x2). For Si(001) 
surface, the p(2x2) structure observed for positive Vb 
is actually attributed to a similar electrostatic effect on 
the basis of density functional theory (DFT). 18 However, 
the observed different dependence on Vb between the two 
types of the kink in Figs. l(a,b) cannot be explained us- 
ing only electric field. Moreover, the calculations by DFT 
under electric field are controversial. 19 Another possible 



origin is the atomic force from the STM tip. A long-range 
attractive force between the tip and the surface 20 causes 
a local lattice deformation of the Ge substrate over a 
few tens nm in diameter, and can modify the potential 
barrier of the dimer flipping. 

In conclusion, we demonstrate that the hot carriers 
injected from the STM tip to the surface states induce 
systematic motions of surface dimer atoms apart from 
the tunneling point. In the case of the electron injection, 
the critical distance of the kink motion is anisotropic, 
and exceeds 80 nm in the dimer-row direction. It also de- 
pends on whether the surface superstructure is p(2x2) or 
c(4x2). These measurements illustrate the dynamics of 
the injected hot carriers in the surface electronic states. 
The observed anisotropy is consistent with the surface 
electronic structure. The semi-classical model of the hot 
electron propagation in the dimer-axis direction can re- 
produce the observed L c which linearly depends on the 
width of the c(4x2) area between the B-type kink and 
the tunneling point. 
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